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Abstract. The main idea and procedure are presented 
for a new synthetical method on constraining 3-D struc- 
ture of emission regions of pulsars. With this method 
the emission regions can be synthetically constrained by 
fitting multi-wavelength observations features, e.g. pulse 
widths, phase offsets between radio pulse profiles and high 
energy light curves, and radio polarization properties. The 
main technique is based on numerically calculating the 
emission directions along any open field line via including 
aberration and retardation effects under a certain mag- 
netic field configuration, e.g. vacuum dipole field. It can 
be easily extended by involving some further effects or 
magnetic field line configuration. The role of observed lin- 
ear polarization properties in the method is discussed. The 
future application to known X/gamma-ray and radio pul- 
sars will be helpful to test and improve current theoretical 
models on pulsar emission. 



1. Introduction 

Although pulsars have been studied for 40 years, the struc- 
tures of radio and gamma-ray emission regions of pul- 
sars remain controversial. It is known that radio emission 
heights are usually less than 10% light cylinder radius for 
normal pulsars (e.g. Cordes 1978, Ranin 1990, Blaskiewicz 
et al. 1991, Kijak & Gil 1997), but the trans-field struc- 
ture of radio emission regions, both in polar or azimuthal 
directions, were not well identified from observations in 
literature. 

Finding out the structure of emission regions is helpful 
to the study on some theoretical issues, e.g. the structure 
of radio emission beam and origin of gamma-ray emis- 
sion. To explain various kinds of radio average pulse pro- 
files, some authors proposed that radio beam is composed 
of conal and core sub-beams (e.g. Rankin 1983, Gil & 
Krawczyk 1996), while some others suggested that radio 
beam is just composed of randomly distributed sources 
(Manchester 1995, Han & Manchester 2001). Since emis- 
sion beam is formed by radiation of particles from flux 
tubes in emission region, the structure of beam is related 
to the structure of emission region. Therefore, it is not 



clear yet whether the radio emission region is composed 
of conal and core-like flux tubes or randomly distributed 
flux tubes. The controversy for the origin of gamma-rays 
remains in different models, i.e. the slot gap model (e.g. 
Muslimov & Harding 2004), the outer gap model (Cheng 
et al. 2000, Hirotani & Shibata 1999) and the annular 
gap model (Qiao et al. 2004). These models suggest differ- 
ent kind of particle accelerators and different gamma-ray 
emission regions. Therefore, determining the 3-D struc- 
ture of emission regions from observation is necessary to 
test and improve pulsar emission models. 

We are developing a synthetical method to constrain 
the emission regions from observational properties. The 
main technique is numerical calculation for the emission 
directions along any open field line by including aberra- 
tion and retardation effects under a certain magnetic field 
configuration, e.g. vacuum dipole field. In this paper we 
focus on its main idea and some remaining problems. 



2. The new synthetical Method 

2. 1 . The main idea and procedure 

The method uses three basic assumptions. (1) In the co- 
rotation reference frame, the emission direction is assumed 
to be along the tangent of magnetic field line. The I/7 
solid angle of emission cone of a relativistic particle is ne- 
glected for the purpose of simplifying calculation. (2) The 
emission patterns are the same at two opposite poles. (3) 
Vacuum rotating dipolar field is adopted. In the Cartesian 
coordinate system of which the rotating axis is z-axis and 
x-axis locates on the meridian plane, the Cartesian com- 
ponents of the magnetic field vector at any point P(x, y, z) 
are as follows (Cheng et al. 2000) 
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B x = \3xz cos a + sin a ([(3x 2 — r 2 ) + 3xyr n 

+ (r 2 - x 2 )r 2 ] cos[(i? - r)/r c ] + [3x2/ - (3x 2 - r 2 )r n 
-x2/r 2 ] sin[(i?-r)/r c ])}, 
By = {3yz cos a + sin a ( [3^2/ + (32/ 2 — r 2 )r n 

-x2/r 2 ] cos[(i? - r)/r c ] + [(32/ 2 - r 2 ) - 3x2/r n 
+ (r 2 -2/ 2 )r 2 ] sin[(ii-r)/r c ])}, 
£> z = -^r { (3z 2 — r 2 ) cos a + sin a ( [3x2: + 3yzr n — xzr\ 

x cos[(i? — r)/r c ] + [32/2; — 3x2;r n — 2/^n] 
x sin[(i? - r)/r c ])} . 



(1) 



where /i is the magnetic moment of the pulsar, R is stellar 
radius, r is the emission radius from stellar center to P, 
r n = r/r c , r c is the radius of light cylinder, a is inclination 
angle between rotation and magnetic axes. 

The method runs in two steps as shown schematically 
in Fig. 1. Step 1, determine geometrical parameters of the 
pulsar and find out all emission points that can be viewed 
by our line of sight (hereafter LOS). To perform the cal- 
culation, the whole open volume is divided into a number 
of layers from the last open field lines to the magnetic axis 
(see Fig. 1 for the shape of foot points of the field lines in 
each layer on polar cap). (1) Determine or constrain the 
values of a and f, where £ is the viewing angle between 
line of sight and rotation axis. Definite values of a and £ 
will greatly reduce calculation. For simplification one can 
use the values that are obtained by fitting the observed 
position angle (PA) of linear polarization with rotation 
vector model (hereafter RVM, Radhakrishnan & Cooke 
1969). But strictly to say they should be figured out con- 
sistently with the method, because aberration, retardation 
and non-static dipole field configuration all influence po- 
sition angle. (2) Calculate tangents of field lines and then 
modify emission direction by involving aberration effect. 
To do the modification the initial emission direction (tan- 
gential to field line) in the co-rotation reference frame is 
converted to the direction in laboratory reference frame 
via Lorentz transform. (3) Select all the emission points 
where the emission direction is aligned with our LOS. 

Step 2, find out the emission points that contribute to 
emission by fitting the observed multi-wavelength proper- 
ties. The way to do this is based on $ — <j> diagram, in 
which the $ — curves for different layers of field lines 
are plotted, where $ is the pulse longitude of the emission 
from a point on a field line and <fi is azimuthal angle of the 
foot point of the field line. After figuring out the diagram, 
one can place the observed multi-wavelength pulse win- 
dows on the diagram, and then finds out which field lines 
may contribute to the pulse (see lower panel in Fig.l). 
In some cases linear polarization percentage may provides 
further constraints on the emission region. The role of po- 



Step 1 

Find out all possible emission points at open field lines that line of sight can see. 



Calculate tangents of a field line. 
(Emission directions in 
co-rotation reference frame) 



Modification by aberration. 
(Emission directions in 
observer's reference frame) 




Foot points of 
different layers 
of open field 
lines on polar 
cap surface 



Repeat the process for all open field lines 
& find all emission points. 



Step 2 

Select the emission points by reproducing the multi-wavelength observational properties 



Modification with retardation effect. 
Calculate pulse phases for all points. 




Plot - <|) and 4> -r diagrams 
for all possible emission points 
(3>: pulse phase) 
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profiles 



Constrain with the following 
observational properties: 

— pulse widths 

— phase offsets 

— linear polarization percentage 

— polarization position angle 



) (deg.) 



Output the geometrical parameters of emission regions seen by line of sight 

— r : height (from emission point to stellar center) 

— <]> : azimuth angle of foot points of open field lines 

— rj: factor indicating how far the open field line is close to the magnetic axis 



Fig. 1. Illustration to the main procedure of the synthet- 
ical method. 



larization properties in this method is discussed in next 
section. 

For a typical pulsar that has radio and high energy 
emission the observed properties used include (a) the ra- 
dio and high-energy pulse widths, (b) the phase separation 
between radio main pulse and inter-pulse (in the case that 
inter-pulse is observed), (c) the phase offset between radio 
and high energy pulses, (d) radio polarization properties, 
viz. linear polarization percentage and PA sweep, (e) view- 
ing angle determined from the elliptical shape of torus or 
rings in nebula for some pulsars, e.g the Crab and Vela 
pulsars. 

2.2. Comparison with conventional methods 

Comparing with conventional geometrical and relativistic 
models to constrain emission altitudes, our new method 
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has three advantages. (1) It can apply to high emission 
altitudes and complex magnetic field configuration. (2) It 
can reveal asymmetrical structure of emission regions. (3) 
It is more synthetical. Details are as follows. 

Geometrical methods (Rankin 1990, Kijak & Gil 1997, 
Wu et al. 2002) and relativistic methods (Blaskiewicz et 
al. 1991, Dyks et al. 2004, Gangadhara 2005) are two 
main conventional approaches to constrain the radio emis- 
sion heights. In geometrical methods geometry relations 
of emission beam and pure dipole field are used, while 
in relativistic model the phase-delay-to-radius relation is 
used by considering the aberration and retardation effects 
in static dipolar field. But geometrical relations and ap- 
proximations of aberration and retardation effects are only 
valid when emission heights are not too large. However, in 
numerical method such limitation vanishes, because re- 
duction of aberration and retardation effects in complex 
field line configuration can be easily done by vector oper- 
ations. Therefore, numerical method can be used to con- 
strain the gamma-ray emission region and radio emission 
region of some young pulsars where the emission height 
may be close to light cylinder. 

A common feature of conventional geometrical and rel- 
ativistic methods is that the emission region is assumed 
to be symmetrical with respect to the meridian that con- 
tains the rotation and magnetic axes. However, the effects 
of aberration, retardation and sweeping back of magnetic 
fields probably make the emission region seen by LOS 
asymmetrical to the meridian. In our method the assump- 
tions are relaxed, so that asymmetrical properties of emis- 
sion regions may be revealed. 

From intuition, more information may be revealed by 
constraining emission regions synthetically with multi- 
wavelength emission properties than with only radio or 
gamma-ray properties. It has been shown that some in- 
teresting results can be obtained for radio and gamma-ray 
emission regions of PSR B1055-52 by synthetically repro- 
ducing the radio and gamma-ray properties, even with a 
modified geometrical method (Wang et al. 2006). The new 
synthetical method is hopeful to go further. 

3. Discussion 

Besides the pulse width and phase offsets, linear polar- 
ization is another key factor in this synthetical method. 
What results one finally gets depend on how to use the 
polarization data. The polarization can be used in two 
ways. First, the PA data should be self-consistently fitted 
with the method to solve or constrain a, £ and the range 
of layers of open field lines. For some pulsars with nebu- 
lae, e.g. the Crab and Vela pulsars, the viewing angle can 
be determined from the torus and ring-like structure near 
pulsar. In this case the uncertainty in geometry is much 
reduced. Second, linear polarization percentage (hereafter 
LPP) can be used to constrain the polar and azimuthal 
scales of emission region. (1) In the cases of nearly 100% 



LPP, the emission region should be limited within a thin 
layer of flux tube, otherwise the superposition of radio 
waves from different layers, which have different PAs, may 
cause significant depolarization. (2) In the case that LPP 
is much less than 100%, it can not be simply concluded 
that the emission region is thick in polar scale, because 
depolarization may be caused by either superposition of 
emissions with scattered PAs or superposition of orthogo- 
nal modes (Stinebring et al. 1984). Detailed studies on PA 
distribution of single pulses may be helpful to distinguish 
them. 

A realistic magnetic field is different from the vacuum 
rotating dipole field. Both magnetospheric current and 
general relativistic effect can influence the field configu- 
ration. It would be interesting to test how the constrained 
results depends on field configuration. 

Our future work includes constraining (1) gamma-ray, 
X-ray and radio emission regions of gamma-ray pulsars, 

(2) X-ray and radio emission regions for X-ray pulsars, 

(3) radio emission regions for radio pulsars with high lin- 
ear polarized percentage and good-quality data of position 
angle. 
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